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We explicitly calculate the Green functions describing quantum changes of topology in Friedman- 
Lemaitre-Robertson- Walker Universes whose spacelike sections are compact but endowed with dis- 
tinct topologies. The calculations are performed using the long wavelength approximation at second 
order in the gradient expansion. We argue that complex metrics are necessary in order to obtain 
a non-vanishing Green functions and interpret this fact as demonstrating that a quantum topology 
change can be viewed as a quantum tunneling effect. We demonstrate that quantum topological 
transitions between curved hypersurfaces are allowed whereas no transition to or from a flat section 
is possible, establishing thus a selection rule. We also show that the quantum topology changes 
in the direction of negatively curved hypersurfaces are strongly enhanced as time goes on, while 
transitions in the opposite direction are suppressed 
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I. INTRODUCTION 

Combining the principles of General Relativity with 
those of Quantum Mechanics seems necessary in order to 
understand more completely the physical nature of the 
gravitational field. As a result, new properties are ex- 
pected to emerge. For instance, in classical General Rel- 
ativity, topology changes are forbidden in the sense that 
their presence would necessarily imply the appearance 
of either singularities or closed timelike curves, a result 
known as the Geroch theorem, see Ref. On the other 
hand, topology changes are widely believed to become 
allowed if the quantum-mechanical nature of the gravita- 
tional field is taken into account. Studying this problem 
is notoriously known as a difficult technical task and in 
the literature only very general aspects of the issue have 
been discussed so far, see e.g. Refs. 0,0. Recently, a 
model where the theoretical ideas about topology change 
can be implemented concretely, at the level of equations, 
has been proposed in Ref. 0. The metric of this model 
is given by a Friedmann-Lemaitre-Robertson- Walker-like 
metric where the curvature of the spacelike sections, usu- 
ally denoted as /C, is now a time dependent function, 
K = K(t) 



ds z = -N\t)dr + a\t)\ d£ 



sin(yA^) 



(1) 

The three-dimensional spacelike hypersurfaces V 3 of the 
model may be endowed with a large class of topologies 
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compatible with such geometry. Technically, the metric 
has always the form displayed above (since it is locally 
defined) but the ranges of variations of the coordinates 
are modified. The scalar curvature of the metric is al- 
ways constant (for fixed t) and equal to 6/C/a 2 . We will 
restrict our considerations to compact (and orientable) 
spaces in order to avoid possible surface terms in the 
Hamiltonian formalism. We can show that any compact 
three manifold with constant curvature is homeomorphic 
to V 3 /r, where the universal covering V 3 is either M 3 , 
S 3 or H 3 , according to the sign of K (K, = 0, K < 0, 
K. > 0, respectively). The group T is the group of cover- 
ing transformations. In three dimensions the closed (i.e. 
compact without boundary) space is now homeomorphic 
to a polyhedron the faces of which are identified by pairs, 
for a general review see Ref. 

Let us first consider the case K. > 0. The universal cov- 
ering is § 3 . Since it is compact, all the three-dimensional 
spaces admitting this universal covering are also com- 
pact. There is an infinite number of such spaces. We 
will only consider two cases: S 3 itself and the Poincare 
dodecahedral space, D 3 = § 3 //*, where I* is the binary 
symmetry group of the icosahedron. The ranges of vari- 
ation of the coordinates (£,,6,tp) can always be written 
as 



< £ < 



S(0,^;V 3 ) 



< < 7T < If < IT. (2) 



The symbol V 3 in the argument of the function 5 indi- 
cates that S(#, ip; V 3 ) is not the same function for differ- 



ent 



(S 3 or 



Of course, the case S is very simple 
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and we just have S(#, ip; § 3 ) = tt. This guarantees that no 
conical singularities will appear in this case. If V 3 = D 3 
the function S is also not arbitrary but its explicit ex- 
pression is much more complicated. For the case K, < 0, 
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where the universal covering is H 3 , the classification of 
closed three-dimensional spaces is still an open question, 
but it is known that the volumes of the polyhedron are 
fixed as in the K. > case. Finally, the case K, = 
have six different tessellations by polyhedrons, some of 
them with arbitrary volume for their polyhedrons (for 
more details see again Ref. What is important for 
us here is that the topologies compatible with one /C are 
not compatible with other K. of distinct sign. Hence, a 
change in the sign of or in the sign of the curvature 
scalar 6/C/a 2 , necessarily indicates a change of topology. 
Our strategy here is to investigate whether this change 
is possible quantum mechanically. 

The quantization of this model was carried out in 
Ref. . In that article, it was shown that a Hamiltonian 
formulation necessarily requires a passage to a midisuper- 
space description. It is not possible to construct a min- 
isuperspace Hamiltonian from the metric Eq. JQ) because, 
as far as K. depends on time, this metric does not rep- 
resent a spatially homogeneous space-time, in the sense 
that the components of the four dimensional curvature 
tensor in a local frame are functions of t and £. The 
existence of the non-vanishing component of the Ricci 
tensor, R t ^ = {£k(t)]/[a(t)N(t)] ^ 0, is a consequence of 
this fact. Hence, we were forced to introduce a midis- 
uperspace model having a non-vanishing shift function 
N^(^,t). The metric was written as 



ds 2 = 



-iV 2 (£,i) + 



N^,t) 



At 2 + 2iV e (£,*)d£dt 



-a 2 (Z,t) 



(3) 



The first step was to carry out the quantization of this 
midisuperspace model. Then, in a second step, we took 
into account in the quantum solutions the restrictions on 
the variables a and a which, from Eq. JSJ, allows us to 
recover the metric Q . Consistency requires that we first 
treat the midisuperspace problem in order to come back 
to the minisuper space model afterward. Correspondingly, 
the Wheeler-De Witt equation of the model remains a 
functional differential equation rendering the findings of 
general exact solutions problematic. In Ref. only two 
particular semi-classical solutions were found for which 
the possibility of a topology change at the quantum level 
was explicitly demonstrated. 

The present article aims at investigating the quantum 
behavior of the system described by the metric (JTJ by 
generalizing the midisuperspace description to a full su- 
perspace description in order to get more general conclu- 
sions than those reached in Ref. |4| . 

Quantities of great interest for this purpose are the 
Green functions since they completely characterize the 
quantum evolution of the system. For a point particle, 
the Green function G(xf,tf, Xi,ti) represents the prob- 
ability amplitude to find the particle at point Xf and 
time tf, knowing that it was at point Xi at the initial 
time i;. In fact, it is necessary to smooth out the Green 



function because, defined as before, it is not square in- 
tcgrablc. Physically, this is due to the fact that the 
particle can never be localized exactly. In the present 
context, a Green function can also be defined and a sim- 
ilar interpretation can be made. In that case, the wave 
function depends on the volume of the spacelike sections 
y (x a 3 , on the three curvature R, and on a dust field 
x{t) (rigorous definitions are given below), which plays 
the role of time: * = ^(y,R,x), see Ref. @- The Green 
function G(%jf,Rf, Xi'iUh^-UXi) can now be defined as the 
probability amplitude of having a volume yi and a three- 
curvature Rf at "time" Xf knowing that, at initial Xh the 
space volume was y\ and the three-curvature R\. There- 
fore, the Green function defined above represents the am- 
plitude of probability for a topology transition as soon as 
Rf has not the same sign as R[. This is the quantity that 
we would like to calculate. 

However, a full calculation of the Green functions is 
probably beyond our present computational capability 
and therefore, we are forced to rely on some approxi- 
mation methods. Firstly, we will restrict ourselves to 
a semi-classical approximation of the Wheeler-De Witt 
equation. In this framework, the wave function will be 
given by the semi-classical wave function ^ — exp(iS/h). 
In the case of a topology change, the phase S has 
a non- vanishing imaginary part, as in ordinary quan- 
tum tunneling, because topology change cannot be ob- 
tained classically, see Ref. 0. Secondly, we will as- 
sume that the spatial sections of our Universe model 
are made of compact homogeneous backgrounds together 
with small perturbations around them. As we are in- 
terested in large scale quantum changes of topology 
in the homogeneous background and not on the small 
scale ones, typical of Wheeler's space-time foam, we 
will use the l ong - wavelength approximation, developed 
m Refs. HH 10, 11], which is valid whenever the char- 
acteristic scale of the problem is much bigger then the 
Hubble radius, and permits a complete analysis of the 
semi-classical wave function and Green functions in the 
full superspace. This will allows us to exhibit explicit 
analytical expressions for the Green functions. In partic- 
ular, we will demonstrate that, for some transitions, the 
corresponding Green function vanishes, establishing thus 
selection rules for topology changes. 

This article is organized as follows. In the second sec- 
tion, we briefly review the long wavelength approxima- 
tion and the evaluation of the Hamilton-Jacobi function 
up to the second order (first order in the curvature). 
Then, in the third section, we perform the calculations in 
the case of real metrics and show that there is no quan- 
tum change of topology in this case. We also present 
a generalization of the results obtained in Ref. [l(| • In 
the fourth section, we argue that complex metrics are 
necessary in order to obtain quantum changes of topol- 
ogy. We explicitly calculate the corresponding Green 
functions and discuss the main consequences that can 
be drawn from their expressions. In the last and fifth 
section we present our conclusions. 
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II. THE LONG WAVELENGTH 
APPROXIMATION FOR QUANTUM GRAVITY 
WITH DUST 

In this section, we discuss the long- wavelength ap- 
proximation technique used to solve the Hamiltonian- 
Jacobi equation which is nothing but the equation for the 
phase functional in a Wentzel-Kramer-Brillouin (WKB) 
semi-classical approximation to canonical quantum grav- 
ity. This technique consists in the expansion of the phase 
functional (or the generating functional in a Hamiltonian 
formulation of General Relativity) in a series of spatial 
gradients and was introduced by Salopek, Stewart and 
Parry 0,0, EJ as a method for solving the full Hamilton- 
Jacobi equation for gravity interacting with matter. In 
Ref. [lfj these authors were able to exhibit a recurrence 
formula from which they could obtain solutions to the 
Hamilton- Jacobi equation to any order of approximation 
for matter made of scalar and dust fields. This will be 
crucial for our treatment because we shall use the second 
order solution to obtain the Green functions associated 
with the WKB wave functional. As we will show, the so- 
lution in the second order approximation will allow us to 
describe changes of topology because the terms in the 
second order expansion of the WKB phase functional 
explicitly contain spatial curvature terms. As already 
mentioned, changes in the sign of these terms indicate a 
change of topology. 

Let us now briefly describe the formalism introduced 
in Refs. [H, 0, The action for General Relativity 

interacting with a dust field x is 



S = 



x V-.9 



2 k 



-Mr- 



2m 



-v(x) 



(4) 



where k = 87r/m pl , mpi being the Planck mass. ^R is 
the Ricci scalar of the space-time metric g^, and \ is a 
velocity potential for irrotational dust particles with rest 
mass m. V(x) is a potential for the dust field x an d 71 is 
the rest number density of the dust particles. The dust 
field x defines the four velocity field for dust particles by 



m 



(5) 



so that x — const, will determine a congruence of space- 
like hypersurfaces foliating the space-time (notice that 

'X,ij.X,v — —m 2 ). Therefore, the dust field can be 
used as the time variable for our model, as well as in the 
Schrodinger-type equation obtained from the Wheeler-De 
Witt equation. 

In the Arnowitt-Deser-Misner formalism the line ele- 
ment reads 



ds 2 



{-N 2 + 1 i3 N l N j )dt 2 +2N i dtdx l +~ flj dx i dx l , (6) 



where N and iVj are the lapse and shift functions, re- 
spectively, and 7ij is the three-metric of the spacelike 



hypersurface. Then, the action (0J can be rewritten as 
S = [ d 4 x(^ X + n ijA fij -NH- N'Hi) , (7) 



where 7r y are the momenta conjugate to 7^ and -k x — 
n r y 1 / 2 {l + x.iX'V™ 2 ) 1 ^ 2 is the momentum conjugate to 
the dust field x- The quantities H, Hi are given by the 
following expressions 

- ±-l 1/2 R + (m 2 + x,iX l Y /2 * x + l 1/2 V{x) ,(8) 



Hi 



(9) 



where R is the Ricci scalar of the three-metric 7^ . 

At this point it is well worth making some dimensional 
analysis, in order to connect the conditions of validity of 
the long-wavelength approximation to the Hubble radius 
of the Universe. From Eq. J5J), we can see that the long- 
wavelength approximation implies that the first term on 
the right-hand side is much larger than the second term, 
at first order in the approximation. From the expression 
for TT ij = (K ij - h i3 K)/K, see Ref. 0, we can see that 
7r^ has the dimension L~ 3 (in the system of units where 
h = c = 1) where L stands for a length. More explicitly, 

1 



4lV 



where £p\ — 1/ mpi is the Planck length and £-r is essen- 
tially the Hubble radius. Therefore, the validity of the 
long- wavelength approximation imposes that, see Eq. JHJ, 

where the wavelength A is a typical length associated with 
the scales where spatial gradients are relevant. In other 
words, A must be greater than the characteristic Hubble 
radius of the model. In our scheme, topology change is a 
large-scale (global) phenomenon and the long-wavelength 
approximation is therefore the appropriate approxima- 
tion to probe it. The above condition, obviously, should 
not contain the Planck length because Eq. JSJ is a classi- 
cal equation. Furthermore, in principle, there is no prob- 
lem in applying this condition to closed models, even 
when their volumes becomes close to the Planck volume, 
(where, anyway, our semi-classical approximation should 
not be valid) if the Hubble radius at this epoch is smaller 
than the Planck length. 

At the quantum level the above system may be quan- 
tized following the Dirac's prescription ^^l- The super- 
Hamiltonian and super-momentum constraints given by 
Eqs. JSJl and © become operators and, when applied to 
the wave functional of the system, result in two relations 
which express that only a restricted region of the Hilbert 
space of wave functionals contains the physical states of 
the theory: 



(10) 
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The first equation is the well-known Wheeler-De Witt 
equation whereas the second one is the so-called quantum 
momentum constraint. In what follows, our treatment of 
quantum changes of topology will be done at the level of 
the semi-classical approximation. This means that the 
wave functionals are assumed to be of the WKB form, 



namely 

<I> = e iS/h , (11) 

where S is the action. At order K°, Eqs. JSJ and @ 
reduce to the Hamilton-Jacobi functional equation 



H(x) = j~ 1/2 K 



SS 



ss 



[2ju (%hjk (x) - Jij (x)^f M (x)\ + \Jm 2 + Y j X,iX,j 



r 



SS 

Sx(x) 



r 1/2 vu)-f- 7 1/2 i? = o, 

(12) 



and to the momentum constraint equation 



Hi(x) 



SS 



S"fkj(a 



SS 



S"fnk(x) 



Ilk,' 



ss_ 

'6X~ X ' 3 



= 



(13) 

the solution of which is the phase S of the WKB wave 
functional, see Eq. JT]J. In Ref. 0, Parry, Salopek 
and Stewart have been able to derive solutions for these 
equations using the long wavelength approximation. The 
method is based on the so-called spatial gradient expan- 
sion of the phase functional S. It consists in expanding 
S in a series of terms according to the number of spatial 

gradients they contain, namely S = S ,(0) +S' (2) +5 (4) H . 

As a consequence, the Hamilton-Jacobi equation can be 
grouped in terms with an even number of spatial deriva- 
tives H = H (0) + H (2) + H (4) + • • •. At zeroth order, the 
Hamilton-Jacobi equation reduces to 

7 *'"« — — — ( 2 7ifc7& - lijlM) 



n (o) 



+m^— + 7 1/ Mx) = 0. 
OX 



(14) 



For a dust field x with a vanishing potential, V(x) = 0, 
the invariance under diffeomorphisms of the generating 
functional suggests a solution of the form 

2 f d 3 x^ 2 H(x). (15) 



This functional satisfies the Hamilton-Jacobi relation if 
H(x) is a solution of H 2 = ~(2m/3)dH/dx yielding 

2m 



(16) 



3(x - X) ' 

where x is an integration constant. The functional S^ ' is 
clearly invariant under coordinate transformations and, 
therefore, satisfies Eq. (|13fl . In addition, it contains no 
spatial derivatives. For the class of generating function- 
als given by Eq. I|15fl . the second-order Hamilton-Jacobi 
equation reads 

SSM 1 H 6S<® 

+ — l 3 X,iX,j 



H 



(2) 



2H 1k 
1 

~2k 



5S& 
Hi 

7 1/2 i? = 



2m 







(17) 



In Ref. |T3, Parry, Salopek and Stewart obtained a dif- 
feomorphisms invariant solution for S^ 2 ' which can be 
written as 



S (2) = I 
K 



d 3 x 1 1 / 2 J( x )R 



(18) 



where the function J(x) is obtained by substituting 
Eq. (IBJ} into Eq. and 

using the zeroth order solution 

(HHJ). This leads to 



dJ 

HJ + m— = - 



dx 



1 

2' 



which yields for S^ 2 ' 



S^ = - 

K 



d 3 x 7 1/2 



10m 



(x - x) + D(x - x) 



(19) 



R 



where D is an arbitrary integration constant. We no- 
tice the fact that S^ contains the three-curvature term 
R whereas this is not the case for S^°\ Since a topol- 
ogy change is signaled by a change of sign of the three- 
curvature R, it is necessary to push the long wavelength 
expansion to the second order in order to have non trivial 
effects. In other words, the quantum topology changes 
can only be studied at this (second) order. 



III. THE GREEN FUNCTION WITH 
LORENTZIAN METRICS 

A. General Derivation 

Let us now consider the semi-classical wave function up 
to the second order in the long wavelength approximation 
* = exp{i[S^ + S<&]/h} with and given by 
Eqs. (|15l) and l|2U|) . as explained in the previous section. 
We restrict ourselves to geometries of the form given in 
Eq. Q. Then, the integrals in Eqs. (JTHJ and 1(20)1 can be 
performed exactly yielding 



^(y,R,x) = ex p|^ 

yRD 



-2yH 



3yR{x - x) 

10m 



(x-x) 2/3 



(21) 
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where the quantity y given by 

dW /2 = Va 3 



y = a x-y 



(22) 



is the total volume associated with the metric (JIJ. V 
is the constant volume of the cell corresponding to the 
homogeneous polyhedron of the tessellation of the spatial 
sections that are compactified. As previously, H can be 
expressed as H = 2m/[3(x — X.)] an d R is the three- 
curvature which can be written as R = 6/C/a 2 . The 
Green functions G are defined by 

/>0O POO 

*(y,R, X ) = / dy' dR'G(y, R, X ; y, R', x) 

JO J -oo 

x*(y',R', X ')- (23) 

Our goal is now to establish the expression for the Green 
function itself from the above equation. For this purpose, 
it is convenient to use a new variable v, which is a combi- 
nation of y and R, namely v = yR. To obtain the Green 
function, we use the fact that Eq. ((23(1 is valid for any 
value of the integration constants x an( i D. However, 
instead of the integration constants x an( i D, it will be 



more appropriate to work with new integration constants 
defined by the following relations (for fixed %') 



r/~\ 



H'{x) 



2m 



Kx'-x) 



(24) 



J'(X,D) ee ^-( X > - x) + D( X ' - x)- 2/3 ■ (25) 



We have just replaced the two constants (x, D) by two 
new constants (H',J'). Let us notice that J' is in fact 
J{x')-> see Eqs. H18fl and l|20|l . hence its name. If we 
inverse the above relations, one arrives at 



X(H') = x' 



2m 

Jit 



, / 2m 



2/3 



3 /2m 
10m \ZH' 



5/3 



(26) 



(27) 



Let us stress again that these expressions are defined for 
a fixed value of x' ■ With these new definitions and the 
explicit expression of the wave function given in Eq. 121fl . 
Eq. l(2l^l now reads 



exp 



- 



-2yH 



3v[x-x(H')} , vD{H',J') 



10m 



\K - X(H')} 2/3 



o y 



x exp 



Hk 



dv'Gl y, -,x;y', —,x' 
y y / 



(—2y'H' + v'j') 



(28) 



We see the effect of having introduced the "new" integra- 
tion constants: the left-hand-side of the previous equa- 
tion now depends on both x an d x' ■ I n order to extract 
G(y,v/y,X',y',v' /y',x')i we multiply the previous equa- 
tion by the factor 



exp 



v J ) 



(29) 



and integrate it over J' and H' from — oo to oo. The re- 
sult of this integration on the right-hand-side of Eq. Ij28() 
is very simple, yielding two Dirac delta functions whose 
arguments are v' — v" and y' — y". This allows us to 
perform the integration over v' and y' immediately. The 
result reads 



G(y,-,x;y", v ^,x') 
\ y y" J 



y" 



27T 2 H 2 K 2 



AH' I dJ'exp<J — 



H' + b 



5H' I V H' + b 



2/3 



V (1 

5 U 



1 

H 1 



(30) 



where 6 = 2m/[3(x — x')]- The first term in the argument H = 2m/[3(x — x)] and then Eqs. 1)26(1 and 1(27(1 . The 
of the exponential has been obtained by using the relation two last terms have been obtained by using the explicit 
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expressions of \ an d D, see Eqs. (|26|) and (|27(l . The above 
equation is the main result of this section. Clearly, the 
central issue is now the calculation of the integrals over J' 
and H'. From Eq. J3UJ), we see that the integral over H' is 
to be performed along the real axis. However, in the next 
subsection, dealing with the case of flat hypersurfaces, 
the integration over H' will be performed in the complex 
plane, along a closed contour. It is important to notice 
that this procedure is just a technical trick, allowing us to 
simplify the calculation of the integral by the use of the 
Cauchy theorem, but it has no deep physical meaning. 
On the contrary, in section IV, the integration in the 
complex plane has a physical meaning and, as we will 
argue, is in fact necessary. Of course, all the following 
complex contours turn out to be, after having taken the 
appropriate limits, equivalent to the integration along the 
real axis as given in Eq. I|30|) . 

B. The Case of Flat Hypersurfaces 



i i i i | i i i i r 
Im(H') 




FIG. 1: The contour of integration C in the complex plane 
[$t(H'), Ss(-ff')]. The contour chosen in this article is the solid 
line contour. Because of the Cauchy's theorem, C is equivalent 
to the small circle contour centered at H' = —b. 



In this section we treat the case where either the start- 
ing hypersurface or the resulting hypersurface (after the 
topology change) is flat. 

Let us first start with the case R = and R" ^ 0, 
hence v — but v" ^ 0. In this case, the integration 
over J' is very simple and yields the Dirac function S(v"). 
Consequently, in order to obtain a non-vanishing Green 
function, one must take v" — 0, that is to say R" = 0. 
This means that the three-curvature remains the same 
(i.e. zero) and there is obviously no topology change. At 
this point, however, we go further and depict the calcu- 
lation of the Green function despite the fact that there is 
no change of topology and despite the fact that the final 
result is already known, see Ref. [TT], The main reason 
for doing this is that this will exemplify techniques that 
will be used in this paper later on. This calculation is 
therefore a warm-up. Eq. (|3(J[1 now reduces to 



G( y ,o,x;v",o,x') 



y 



2tt 2 H 2 K 2 



AH' exp 



i (2ybW 
Jm\H' + b 



2y"H> 



(31) 



As mentioned before, the integral in the right-hand-side 
of Eq. H31|) can be calculated by extending H' to the 
complex plane with the choice of the integration contour 
represented by the solid line in Fig.^ This contour con- 
tains the pole H' = —b of the exponential function ap- 
pearing in the right-hand-side of Eq. I|31|l . The fact that 
H' = —b, which is a pole of the argument of the exponen- 
tial, is also a pole of the exponential itself can easily be 
seen if one uses the Taylor expansion of the exponential 
function. When the limits p — * +oo and e — > are taken, 
see Fig. ^ the integral over H' in the right-hand-side 
of Eq. (|31J) reduces to an integration on the real axis, 
the contribution coming from the big semicircle in the 
upper-half complex plane going to zero as a consequence 



of the Jordan's lemma. Because of the Cauchy's theo- 
rem, the integration along C is in fact equivalent to the 
integration over a small circle centered around H' = —b, 
see Fig. ^ This was the contour considered in Ref. 0] 
which is therefore equivalent to the contour C used in this 
article. As a consequence, one obtains 



G{y,o,x;y",o, x ') = 



y 



I d-ff'exp 
'c 



2ir 2 h 2 K 2 
i f 2ybH' 

Hk\ W + b 



2y"H' 



(32) 



It has been shown in Ref. that the integration over 
H' along the contour C can then be done explicitly. The 
exact result is expressed in Eq. (4.25) of Ref. [Tj in terms 
of a Bessel function of order one. We do not present this 
result here since this aspect of the calculation has been 
given in details in Ref. . 

Let us now consider the other possibility where R ^ 0, 
hence v ^ but now R" = (v" = 0). Despite the 
fact that the argument of the exponential function is now 
more complicated, the integration over J' can still be 
performed very easily and yields the Dirac delta function 



H' 



2/3' 



(33) 



Since the argument of the Dirac function has no zeros in 
H', the integration over H' leads to a vanishing Green 
function. 

Hence, we have established the following result: there 
is no quantum change of topology from or to a flat space- 
like hypersurface in the homogeneous background, at 
least at second order in the long- wavelength approxima- 
tion. 
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C. The Case of Curved Hypersurfaces 

We now assume that both R and R', and hence v and 
v", do not vanish. This time the integration over J', in 
the right-hand-side of Eq. (|3U|) , gives the following Dirac 
function [to be compared with Eq. 



2/3 



(34) 



Since H' is real, the argument of the Dirac function has 
zeros if and only if R and R 1 (or, equivalently, v and v") 
have the same sign and, in this case, the zeros are given 
by 



H' 



n o 



it 



3/2-1 



(35) 



If R and R' have not the same sign, the corresponding 
Green function is zero and quantum topology change is 
forbidden. Therefore, as in the flat case, there is no quan- 
tum change of topology. 

Despite this fact, and for the sake of completeness, we 
now present the resulting Green function since this repre- 
sents a non trivial generalization of the results obtained 



in Ref. ^lj. Let us notice that this case is different from 
the one given in Eq. (|31|) . There, we had a Dirac function 
of v" only and, therefore, it was possible to put the Dirac 
function outside the integral over H'. Here, we face a dif- 
ferent situation since, this time, the Dirac function does 
depend on H' (the second flat case , i.e. v ^ 0, also pre- 
sented this feature but since the argument of the Dirac 
function had no zeros, this case was in fact trivial). In 
order to calculate the integral l|30|) for the present case, 
we express the Dirac function as 



2/3 



36 

2\v\ 



W + b r 
S(H' - H'o) + 5{H' - H'+ 



-5/2 



(36) 



and the integration becomes trivial. As a consequence, 
the Green function is given by the sum of two contribu- 
tions coming from the two Dirac functions in Eq. (|36|l 



G r = (G+ + G-)e 



K_ 



(37) 



where 9 is the step function and where G± are given by 



G+ 
G- 



3my x 



" T 5/2 



ttHk \ v\(x- x") 

imy"x 5 / 2 
irhn | v | (x - x") 







exp< 








exp< 


Hk 



2yb{x 3 / 2 - 1) 
^3/2 

2yb{x a / 2 + 1) 



x 



3/2 



r 



2y"b(x 3 ' 2 - 1) 
2y"b(x 3/2 + 1) 



^1/3(^5/2 _ 



5b(x 3 / 2 - l)x 



hb{x 3 ' 2 + l)x 



(38) 
(39) 



with x = (y/y") 1 / 3 - This completes our study of the 
Green functions in the case of real metrics. We have 
demonstrated that there is no quantum change of topol- 
ogy in that long as this transition involves 
a change of sign in the curvature of the spatial sec- 
tions. Our formulation does not allow to signal a change 
of topology when the transition involves compactified 
spacelike sections with distinct topology but having the 
same sign of curvature, a fact that would require a com- 
pletely distinct and much more elaborated approach to 
the problem. These conclusions are valid up to second 
order in the gradient expansion. 



IV. THE GREEN FUNCTION WITH COMPLEX 
METRICS 



We have just seen in the previous section that there is 
no change of topology if H' remains real. This is linked 
to the fact that the equation v" = v[b/(H' + 6)] 2 / 3 has no 
solution on the real axis of H' if v and v" have opposite 



signs. It is worth recalling again that the extension of 
H' to the complex plane made in section III, in the case 
of flat hypersurfaces, was a mere mathematical trick to 
calculate the integrals, albeit all the poles inside the con- 
tours were located on the real axis. On the other hand, 
the above equation possesses a solution in the complex 
plane, namely 



H' = H' Q 



b 1 + i 





3/2x 


V 


V 1 





(40) 



and, therefore, this strongly suggests to consider H 1 as a 
true complex variable. As we will discuss in the present 
section, the complexification of H', with its natural con- 
sequence namely a complexification of the metric itself, 
is a physical requirement. This will lead us to inter- 
pret topology change as a quantum tunneling effect. In 
the case of complex metrics, the corresponding Green 
function needs to be modified and its expression can be 
obtained from the expressions derived before. We now 
describe how this can be done. 
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A. General derivation 

Changes of topology described by the metric JIJ re- 
quires a complicated midisuperspace formulation of the 
problem, as described in the Introduction and explained 
in more details in Ref. p|. However, in the case of 
curved spatial hypersurfaces, the midisuperspace formu- 
lation can be circumvented and the problem greatly sim- 
plified by the introduction of complex metrics. Let us 
now define the metric according to 



ds 2 = -N 2 (t)dt 2 + a 2 (t) 



dz 2 + sin 2 z(d6> 2 + sin 2 9d(p 2 ) 



(41) 

and allow N c and z to be complex. Setting N c = N 
and z = £ to be real yields, let us say for the final 
configuration, a positively curved spatial hypersurface 
in a Lorentzian four-geometry, while requiring initially 
that N c = —iN and z = — i£ are purely imaginary 
yields a negatively curved spatial hypersurface, also in 
a Lorentzian four-geometry [of course, one can inter- 
change the transition by choosing sinh instead of sin in 
Eq. H41f> ]. Therefore, the transition from positively to 
negatively curved spatial sections may be obtained with 
a mere complexification of the metric, without changing 
the functional form of the metric in Eq. (|41|l or with- 
out needing to introduce a time dependence on as in 
Eq. (JTJ. In other words, the transition described above 
can be actually realized by following trajectories in the 
space of complex metrics. Notice also that, with the com- 
plexification procedure, there is an overall change of sign 
in the 4-metric while the transition from the initial to the 
final configuration occurs , implying a passage through 
degenerate metrics. This is typical of transitions with a 
change of topology as pointed out in Ref. [l3|. Let us 
also emphasize that the use of complex trajectories in 
order to describe the tunneling effect is standard and has 
been used elsewhere, for instance in the study of chaotic 
potentials or in the calculation of the rotational spectra 
of molecules, see Ref. [l4j. 

The change in the lapse function, N c = —iN, is equiv- 
alent to a Wick rotation to the imaginary time r = —it, if 



one understands t as the proper time of observers whose 
trajectories are orthogonal to the spacelike hypersurfaces. 
In this context, one needs to adapt the ADM formal- 
ism presented at the beginning of this article in order to 
obtain the Green function that will describe a topology 
change. We now turn to this question. 

The Euclidean action defined as I = —iS[g^\, see 
Ref. [l3j, is given by 



I 



d x y/g^ 



1 

2k 



(4) 



2m 



v(x) 



(42) 



where (?e denotes the Euclidean 4-geometry whose line 
element can be expressed as 

ds 2 = (N 2 + ^N i N J )dt 2 + 2N i dtdx i +j ij dx i dx i . (43) 

Note that we have used the relation \J—g c = —i^/ffE, 
where g c is a Lorentzian metric complexified as in 
Eq. (|41|l . Having determined what the action is, we can 
study the Hamiltonian formalism. Calculations very sim- 
ilar to those performed in the standard ADM formalism 
show that the Hamiltonian obtained from the Euclidean 
action given above reads 

H E = NH E + N'H, , (44) 

where TLe and TLi are respectively given by 

He = -k 1 - l l 2 ^n kl {2 1]klu - lijlu ) - ^-j 1/2 R 

Ik 

+i(m 2 + X,jX' j ) 1/2 ir x + 1 1/2 V{ X ) , (45) 
Hi = -2(r fie n tk ) >k + * tk >ytk,i + * x X,i, (46) 

and where, as before, -k ij are the conjugate momenta 
to jij and 7r x is the conjugate momentum to the dust 
field x- The quantity R is the Ricci scalar of the three- 
metric "fij . From the expression of 7Ye derived above, one 
can establish the Euclidean Hamilton- Jacobi equation. It 
reads 



~,- 1I2 k 



5S 



5S 



[2^ u (x)^ jh {x) - (x)j ki (x)] + i\Jm 2 +^x,iXd 



while the momentum constraint equation remains the 
same as in Eq. (j*H?|l . Comparing Eq. (|17J> with Eq. (fT^l , 
one can see that going to the imaginary time (by means 
of the Wick rotation) changes the sign of the kinetic term 
in the gravitational sector. This is a well-known modi- 
fication when going from Lorentzian to Euclidean signa- 
tures. On the other hand, the matter sector, containing 



the term SS/Sx, is multiplied by complex number i. This 
latter change was also expected since \ pl avs the role of 
time. Setting r = — it implies \c = — ix after the Wick 
rotation, hence the above result. 

Having the Euclidean Hamilton-Jacobi equation and 
the momentum equation at our disposal, we now seek 
solutions using the gradient expansion. In the following, 
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quantities with the index c refer to the metric (|41() after 
the Wick rotation (N c = —iN, z complex) while quan- 
tities without index refer to the (real) Lorentzian metric 
(|41|l (iV=real, £ real). As it was the case before, see 
Eq. I|15|l . the invariance by reparametrization suggests a 
solution of the form 



~ J d 3 x 7 y 2 i? c ( X ). 



(48) 



Inserting this form into the Euclidean Hamilton- Jacobi 
equation yields, at the zero order, H 2 = (2im/3)dH c / 'd\- 
This equation should be compared with the equation 
written after Eq. (|15|) . The difference between those 
two equations is a factor — i as expected. The corre- 
sponding solution can be written as H c = —iH, where 
H(x) — 2m/[3(x — x)]- Again, this solution bears some 
close resemblance with the solution of Eq. fTSjl . 

At second order, the diffeomorphism invariant solution 
reads 



SP = i J d 3 x 7c V2 j c[x)Rc 



(49) 



and is similar to the solution of Eq. (|2f)|> . Applying the 
same method as before, i.e. inserting the above equation 
into the Euclidean Hamilton- Jacobi equation, yields the 
following relation 

TT T 9J C TT dJ C 1 , ^ 

— H C J C + im— — = iHJ c + im— — = - . (50) 
dx dx 2 v ; 

Since we have already established the expression of H c , 
the above relation is a differential equation for the quan- 
tity J c only. The solution is very similar to the one found 
in section II of this paper. It reads J c = — iJ, where, as 
can be deduced from Eq. 1|20[) . one has 



10m 



ix- X) + D{x~ X.) 



-2/3 



(51) 



We now have at our disposal the Euclidean action at 
zcroth and second order in full generality, i.e. for any 
metric tensor. As before, we now restrict our consider- 
ations to the metrics of the form l|41jl. This allows us 



J 



to perform the integrals (|48|l and (|49|l exactly. How- 
ever, in the Euclidean case, one has to pay attention to 
the following fact. When the transformation z — — i£ is 
made, the three-metric gets modified by an overall sign 
(besides the change from sin to sinh). This implies that 
7c' /2 = — ij 1 / 2 and, hence, the volume y c can be written 
as 



d 3 a;7, 



1/2 



-ly . 



(52) 



For the same reason, the three-curvature satisfies R c = 
R{lc V ) — —R{~i^) = —R- Combining the previous re- 
sults, one has v c — y c R c — iyR — iv. 

We are now in a position where we can calculate the 
semi-classical wave function in the long wavelength ap- 
proximation at the time x- The result reads 



*fec-Rc,x) = exp 



exp 



i 



(-2y c H c + v c J c 
(2yH + vJ) . 



(53) 
(54) 



This solution should be compared with Eq. (|21|) . One 
notices that there is a difference, namely a change of sign 
in the first term of the exponential. As we are going to 
show below, this simple change of sign turns out to be 
crucial in order to obtain a Green function which allows 
a quantum change of topology. 

The last step is to calculate the Green function ex- 
plicitly. This Green function describes a transition from 
a situation where the wave function is expressed by the 
standard formula, Eq. (|21J) . to a case where the wave 
function is modified by the above-mentioned change of 
sign, see Eq. (|53|l (or vice versa if one wishes to consider 
the opposite transition). As a consequence, in order to 
obtain the Green function, we just need to repeat the 
calculations performed between Eq. ll'Mli and Eq. H32fl . 
bearing in mind that, in the final (or initial) configura- 
tion, the wave function is now described by Eq. I|53|) . We 
obtain the following result 



G i (y,R, X ;y',R',x') = 



y 



AH' 



dJ'exp< — 



-2yH - 2y'H' - v'j' 



jMx ~ X) 

10m 



vD 



(x-x) 2/3 



(55) 



r 



At this point, one should compare the above result with 
Eq. (|3*2*|l . Being given the previous considerations, it does 
not come as a surprise that the sign of the second term in 
the exponential is different. The origin of the difference 
is clearly the modified wave function of Eq. (J53J. Apart 
from this (crucial) difference, the previous expression is 



identical to Eq. I|30|l . Let us also notice that, in order 
to derive the Green function, we have already assumed 
that the contour C reduces to an integration along the 
real axis. This allows us to extract the Green function 
directly. We will show below that this is indeed the case. 
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Im(H') 
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FIG. 2: The contour of integration C in the complex plane 
[$t(H'),Q(H')] in the case of curved hypersurfaces for the 
modified Green function. The contour chosen is the solid line 
contour which contains the complex pole of the Dirac function 
H but avoids the poles of the exponential function, namely 
H' = and H' = -b. 



B. Explicit calculation of the Green function 

Equipped with the modified Green function l|55(l . we 
now return to our main subject, i.e the study of the quan- 



tum topology changes. We now consider the case where 
v and v' have opposite signs, namely sign(w)= — sign(w'). 



The first step is to perform the integration over J'. 
This yields the Dirac <5-function of Eq. i|34[l (except that 
one should replace v" by v'). This function has the root 
given in Eq. lj4*0"l) . i.e. H 1 — H' . The second step is the 
choice of the contour C. In the complex plane, integration 
of a Dirac S- function can be represented by an integration 
along a small circle around the roots of its argument. For 
this reason, the integration over H' will be performed 
along the contour which goes from — p to p on the real 
axis, avoiding the poles of the exponential function, H' = 
—b and H' = 0. The contour is closed by considering a 
large semicircle in the half part of the complex plane with 
negative imaginary part. This contour is represented in 
Fig. |21 It is equivalent to a small closed contour around 
the pole H' = H' Q because of the Cauchy theorem. In 
addition, as p — > oo, the integral along the large lower 
semicircle goes to zero and we are reduced to an integral 
from — oo to +oo, as before. As announced above, this 
last argument justifies the fact that we have extracted 
the Green function from Eq. 155|) . Performing explicitly 
the integration, the Green function for topology change 
reads 



Smy'x 



5/2 



/ 1 

AttHk I v | (x — x') [fin 



4m(y-y>) SJCy^jx 4 - l)( x - x ') 
3(X - X') 10(x 3 + l)x 



exp 



A 1/2 



Bm(yy') 
3(X " X') 



10(x 3 + 1) 



(56) 



where the quantity x has been introduced before. This 
is the main result of this section. 

Let us now check that the above expression possesses 
the correct properties of a Green function. The Wheeler- 
De Witt equation for the homogeneous background is 
given by 



2 k 



(57) 



The Green function Gi must satisfy the semi-classical 
version of Eq. (|5T|) . Recovering the units, Gi has the fol- 
lowing form: Gi = C exp(iS/h) exp (A/ ft). After insert- 
ing this expression in Eq. I|57|l. one obtains an equation 
whose real and imaginary parts, in the limit ft — ► 0, read 



dS 3k 
m— y 

d X 4 y 



os_ 

dy 



dA 
dy 



ICy 1 / 3 
2k 



0(58) 



dA 3k OSdA 
dx 2 dy dy 



Using the expression of Gi derived above, see Eq. JSHJ, 
one can easily check that Eqs. (|58|l and (|59|l are indeed 
satisfied in the long-wavelength approximation. Let us 
notice that terms involving R 2 = (6K./y 2 ^ 3 ) 2 have been 
neglected. Hence, Gi possesses all the good properties 
of a Green function signaling a topology change: it de- 
scribes the transition from a wave function describing an 
homogeneous spatial geometry with positive (negative) 
curvature to another wave function describing another 
homogeneous spatial geometry with negative (positive) 
curvature [see Eq. I|23|) ]. and it satisfies the semi-classical 
Wheeler-De Witt equation in the long-wavelength ap- 
proximation. 



11 



C. Discussion 

In this section, we briefly discuss the main properties 
of the Green function of Eq. (|56|l . 

First of all, it is interesting to notice that, without the 
change of sign obtained in the first term of the exponen- 
tial, see Eq. (|53[) . the Green function obtained through 
Eq. iffify would not satisfy the Wheeler-De Witt equa- 
tion IpTTjl. This change of sign is a direct consequence of 
describing the change of topology by means of complex- 
ification of coordinates, which in turn implies an overall 
change of sign in the four-metric while the transition is 
taking place. This confirms the ideas of Ref. 01 stating 
that, in a topology transition, one should pass through 
degenerate metrics. We notice in passing that Eq. is 
not exactly the classical Hamilton-Jacobi equation for a 
minisuperspace corresponding to the homogeneous met- 
rics considered here. There is an extra term which in- 
duces the quantum change of topology. This is the quan- 
tum potential in the Bohm-de Broglie interpretation lan- 
guage. This extra term comes from the imaginary part A 
of the action and indicates that such a change of topology 
is in fact a quantum tunneling effect. 

Secondly, examining the Green function l|5()|) . one can 
notice that a topology change is very improbable if 
A% = x ~ x' i s very short. As A% increases, quantum 
changes of topology become possible in both directions 
(from positive to negative curvature and vice- versa), but 
when A\ becomes very large, negative to positive cur- 
vature transitions are suppressed while positive to neg- 
ative curvature transitions are enhanced. This suggests 
that negative curvature topologies are preferred. Finally, 
for finite A%, topology transitions between large volume 
spacelike hypersurfaces are very improbable, as expected. 

V. CONCLUSIONS 

We now briefly summarize the main results reached in 
this article. 



We have obtained a Green function describing a topol- 
ogy change by using a double approximation: the semi- 
classical canonical quantum gravity and the long wave- 
length approximation. Our considerations have been re- 
stricted to transitions between homogeneous spacelike 
hypersurfaces with different intrinsic curvatures. In this 
framework, we have been able to demonstrate that tran- 
sitions involving flat hypersurfaces are forbidden. How- 
ever, we have also shown that this is no longer the 
case if the hypersurfaces are curved. In order to obtain 
the corresponding Green functions describing transitions 
between curved hypersurfaces, we have used imaginary 
coordinates (not only timelike coordinate but also one 
spacelike coordinate) in order to circumvent the midisu- 
perspace problem raised in Ref. Q (i.e. the components of 
the four-dimensional curvature tensor are only functions 
of time). The explicit expression of the Green function 
has been obtained by performing an integration in the 
complex plane. The resulting Green function shows that 
topology changes in the direction of negatively curved hy- 
persurfaces are strongly enhanced as time goes on while 
transitions in the opposite direction are suppressed. In a 
finite amount of time, transitions between large spacelike 
hypersurfaces are improbable. 

The formalism developed in this paper can also be ap- 
plied to more complicated models. This will be the sub- 
ject of our future investigations. 
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